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the terms in the perturbative coefficients. It allows us to systematically determine the argument of 
the running coupling order by order in pQCD in a form which can be readily automatized. The new 
method satisfies all of the principles of the renormalization group and eliminates an unnecessary 
source of systematic error. 
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An important goal in high energy physics is to make 
perturbative QCD (pQCD) predictions as precise as pos- 
sible, not only to test QCD itself, but also to expose new 
physics beyond the standard model. In this letter we 
present a systematic method which determines the ar- 
gument of the running coupling order by order in pQCD 
and which can be readily automatized. The resulting pre- 
dictions for physical processes are independent of theo- 
retical conventions such as the choice of renormalization 
scheme and the initial choice of renormalization scale. 
The resulting scales also determine the effective number 
of quark flavors at each order of perturbation theory. The 
method can be applied to processes with multiple physi- 
cal scales and is consistent with QED scale setting in the 
limit iVc — )■ 0. The new method satisfies all of the prin- 
ciples of the renormalization group [1], and it eliminates 
an unnecessary source of systematic error. 

The starting point for our analysis is to introduce 
a generalization of the conventional schemes used in 
dimensional regularization in which a constant —5 is 
subtracted in addition to the standard subtraction 
ln47r — 7b of the MS-scheme. This amounts to redefin- 
ing the renormalization scale by an exponential factor; 
i.e. /x^ = /x?jgexp(5). In particular, the MS-scheme is 
recovered for 5 = ln47r — 7^;. The (5-subtraction defines 
an infinite set of renormalization schemes which we call 



* mojazaOslac. stanford.edu 
t sjbth@slac.stanford.edu 
t wuxg@cqu.edu.cn 



5-72,enormalization {TZs) schemes; since physical results 
cannot depend on the choice of scheme, predictions must 
be independent of 6. Moreover, since all TZs schemes are 
connected by scale-displacements, the /3-function of the 
strong QCD coupling constant a = as /47r is the same in 
any 72.5 -scheme: 

, 00 

The TZg-scheme exposes the general pattern of noncon- 
formal {/3i}-terms, and it reveals a special degeneracy of 
the terms in the perturbative coefficients which allows us 
to resum the perturbative series. The resummed series 
matches the conformal series, which is itself free of any 
scheme and scale ambiguities as well as being free of a 
divergent renormalon series. It is the final expression one 
should use for physical predictions. It also makes it pos- 
sible to setup an algorithm for automatically computing 
the conformal series and setting the effective scales for 
the coupling constant at each perturbative order. 

Consider an observable in pQCD in some scheme which 
we put as the reference scheme TZq (e.g. the MS-scheme, 
which is the conventional scheme used) with the following 
expansion: 
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Po{Q^)=Y.^.iQVfil)a{iioy , (2) 

4 = 

where /ig stands for the initial renormalization scale and 
Q is the kinematic scale of the process. The more gen- 
eral expansion with higher tree-level power in a can be 
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readily derived [2] and will not change our conclusions 
and results. Since results in any TZs are related by scale 
displacements, we can derive a general expression for p 
by using the displacement relation between couplings in 
any T^-^-scheme: 

a(^o) = + E ^ U^.. (-^)" , (3) 

n— 1 ^ ' 

where we used In /ig/Ma = Then p in TZs to order 
reads: 

Ps{Q^) =ro + riai{pi) + (r2 + l3oriSi)a2ip2f 

+ [ri + P2ri5i + 2/3ir2(52 + 3/3or3^3 + 3/3^r2(5^ 

where /zf = Q^e*', the initial scale is for simplicity set to 
Mo = ^'^'i defined ri(f) = r^. An artificial index 
was introduced on each a and 5 to keep track of which 
coupling each 5 term is associated with. The initial scale 
choice is arbitrary and is not the final argument of the 
running coupling; the final scales will be independent of 
the initial renormalization scale. 

In a conformal (or scale-invariant) theory, where 
{ft} = {0}j the 5 dependence vanishes in Eq.(4). There- 
fore by absorbing all {/?,;} dependence into the running 
coupling at each order, we obtain a final result indepen- 
dent of the initial choice of scale and scheme. The use of 
TZs allows us to put this on rigorous grounds. From the 
explicit expression in Eq.(4) it is easy to confirm that 

dps r,, .dps , , 

Scheme-invariance of the physical prediction requires 
that dps/dS — 0. Therefore the scales must be set such 
that effectively /3(a) =0 in the final expression and this 
leads to the conformal series. The scheme-invariance re- 
quirement is a theoretical requirement, which must be 
satisfied at any truncated order of the pertubative series, 
and is different from the physical fact that the all-orders 
expression for an observable is renormalization scale and 



scheme invariant; i.e. dp/dfiQ = 0. This demonstrates 
the concept of the principal of maximal conformality [3] 
(PMC) to any order. 

The expression in Eq.(4) exposes the pattern of {f3i}- 
terms in the coefficients at each order. Such a pattern was 
recently considered in Ref. [4] . The 7?.5-scheme reveals its 
origin and its generality for any pQCD prediction. It is 
possible to infer more from Eq. (4) ; since there is nothing 
special about a particular value of S, we conclude that 
some of the coefficients of the {/3i}-terms are degenerate; 
e.g. the coefficient of /3oa(Q)^ and l3ia{Q)^ can be set 
equal. Thus for any scheme, the expression for p can be 
put to the form: 

P(Q^) =fo,o + rifla{Q) + [r2,o + l3or2,i]a{Q)'^ 

5 

+ [''4,0 + /32?'2,i + 2/3ir34 -I- -/3i/?o''3,2 + 3/3o?'4,i 

+ 3/30^4,2 + Plr^,MQf + Oia^) (6) 

where the are the conformal parts of the perturba- 
tive coefficients; i.e. — r^^o + C({/3i})- The 7?.5-scheme 
not only illuminates the {/?i}-pattern, but it also exposes 
a special degeneracy of coefficients at different orders. 
The degenerate coefficients can themselves be functions 
of {/3i}, hence Eq.(6) is not to be understood as an ex- 
pansion in {/3i}, but a pattern in with degenerate 
coefficients that must be matched. We have checked that 
this degeneracy holds for several known results. 

The expansion in Eq.(4) reveals how the {/3i}-terms 
must be absorbed into the running coupling. The differ- 
ent 6k 's keep track of the power of the I /e divergence of 
the associated diagram at each loop order in the follow- 
ing way; the (5^a"-term indicates the term associated to a 
diagram with l/e"~'° divergence for any p. Grouping the 
different (5fc-terms, one recovers in the TVc — >■ Abelian 
limit [5] the dressed skeleton expansion. Resumming the 
series according to this expansion thus correctly repro- 
duces the QED limit of the observable and matches the 
conformal series with running coupling constants evalu- 
ated at effective scales at each order. 

Using this information from the 5fc-expansion, it can be 
shown that the order a{Q)'' coupling must be resummed 
into the effective coupling a{Qk)'', given by: 
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(7) 



rfe,oa(Qfc)'=^rfc,oa(Q)''+rM fca(Q)''"V(a){^fe,i+Ai'^(a)i?M + --- + Ai" '^(a)i?fc.„} , (8) 
which defines the PMC scales Qk and where we introduced 



oa a 
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Eq.(8) is systematically derived by replacing the In-* Q\IQ'^ by Rk,j in the logarithmic expansion of a{Qk)'' up to the 
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highest known i?fc.„-coefficient in pQCD. The resummation can be performed iteratively using the renormahzation 
group equation for a and leads to the effective scales for an NNNLO prediction^: 

Rk,i+Ai'\a)Rk,2 + A^'\a)Rk.3 ^ ^^^^ 

The final pQCD prediction for p after setting the PMC scales Qi then reads 

p{Q^) ^ro,o + ri,oa(Qi) + r2,oa{Q2f + rs^aiQif + r^^oaiQi)'^ + ©(a^) , (11) 



Note that 04 remains unknown. This last ambiguity re- 
sides only in the highest order coupling constant, and is 
negligible in practice. 

It is easy to see to leading logarithmic order (LLO) that 
the effective scales are independent of the initial renor- 
mahzation scale fiQ. This follows since taking /io =/= Q wc 
must replace Rk^i — -Rfc^i-t-lnQ^Z/iQ and thus the leading 
order effective scales read InQ^ lo/A'o = Rk,i+^^Q^ / fJ'o, 
where /io cancels and Eq. (10) at LLO is recovered. This 
generalizes to any order. In practice, however, since 
the /3-function is not known to all orders, a higher or- 
der residual renormalization-scale dependence will enter 
through the running coupling constant. This residual 
renormalization-scale dependence is strongly suppressed 
in the perturbative regime of the coupling [6] . 

The effective scales contain all the information of the 
non-conformal parts of the initial pQCD expression for 
p in Eq.(6), which is exactly the purpose of the running 
coupling constant. The quotient form of Eq. (10) sums 
up an infinite set of terms related to the known rj^k^Q 
which appear at every higher order due to the special 
degeneracy. It is, however, not the full solution since 
this requires the knowledge of the rj^k^Q-teims to all 
orders. The method systematically sums up all known 
non-conformal terms, in principle to all-orders, but is in 
practice truncated due to the limited knowledge of the 
/3-function. 

In earlier PMC scale setting [3], and its predecessor, 
the Brodsky-Lepage-Mackenzie (BLM) method [7], the 
PMC/BLM scales have been set by using a perturbative 
expansion in a and only approximate conformal series 
have been obtained. Here, we have been able to obtain 
the conformal series exactly due to the revelation of the 
{/3j}-pattern by TZs; the effective scales have naturally 
become functions of the coupling constant through the 
/3-function, in principle, to all orders. 

In many cases the coefficients in a pQCD expression 
for an observable are computed numerically, and the {(3i} 
dependence is not known explicitly. It is, however, easy 
to extract the dependence on the number of quark fla- 
vors Nf , since Nf enters analytically in any loop diagram 
computation. To use the systematic method presented in 



this letter one puts the pQCD expression into the form 
of Eq.(6). Due to the special degeneracy in the coeffi- 
cient of the {/3,i}-terms, the Nf series can be matched to 
the rj^k coefficients in a unique way^. This allows one to 
automate the scale setting process algorithmically. 

The general iV/-series of the n-th order coefficient in 
pQCD reads: 

r„ = Cn,o + c„,iiV/ + • • • + c„,„_iiVJ-i . (12) 

By inspection of Eq.(6) it is seen that there are exactly 
as many unknown coefficients in the {/3i}-expansion at 
the order a" as the Nf coefficients, c„j-. This is real- 
ized due to the special degeneracy found in (6). The rij 
coefficients in Eq.(6) can thus be expressed in terms of 
the Cnj coefficients. This means that the iV^ terms can 
unambiguously be associated to {/3i}-terms and demon- 
strates PMC as the underlying principle of BLM scale 
setting. The relations between c„j- and rij are easy to 
derive and they transform the BLM scales into the cor- 
rect PMC scales [2]. 

The automation process can be outlined as follows: 

1. Choose any (5-7?,enormalization scheme and scale. 

2. Compute the physical observable in pQCD and ex- 
tract the Nf coefficients, cj. j. 

3. Find the /3j coefficients, r^j from the c^.j coeffi- 
cients and compute the PMC scales, Qk- 

4. The final pQCD expression for the observable reads 

Pfinal(Q) = J2k=0'^k,0a{Qk)''- 

Example: e+e^ hadrons. The ratio for electron- 
positron annihilation into hadrons, R'^ was re- 
cently computed to order a* [8] and can be shown to 
exactly match the generic form of Eq.(6). It can be de- 
rived by analytically continuing the Adler function, Z?, 
into the time-like region, with D given by: 

D{Q^)^^{a)^P{a)^^{Q\o) , (13) 



Detailed derivations will be given elsewhere [2] 



^ In principle, one must treat the Nf terms unrelated to renormal- 
ization of the gauge coupling as part of the conformal coefficient. 
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where 7 is the anomalous dimension of the vector field, 
n is the vacuum polarization function and they are given 
by the perturbative expansions: 7(a) — X^i^o'T"'^" 
n(a) = nno". It is easy to show that to order a"^ 

the perturbative expression for W^'' in terms of 7,1 
and n„ reads: 

^e+e-^h(g) ^ ^ ^^^(g) ^ _ /3^U2HQf (14) 

2 

+ [73 + f3iU2 + 2^on3 - (3^^HQf 

+ [74 + /^2n2 + 2/3in3 + 3/3on4 

This expression has exactly the form of Eq.(6) with the 

2 

identification: = 7i: n,! = T^i, n,2 ^ -^lz~2 Siiid 
ri^3 = — 7r^ni_2. The 7^ contain A'^^-terms, but since they 
are independent of 6 to any order, they are kept fixed in 
the scale-setting procedure. Note that we have knowledge 
of even higher order r^ j coefficients, and this allows us 
to set the effective scales Qi, Q2 and Q3 to the NNNLO, 
given by Eq.(lO). It is worth noting that the Adler func- 
tion D itself has a much simpler {/?i}-structure. By con- 
vention the argument of a is space-like; thus the 7r^-terms 
appearing in i?'^^^ could be avoided by using a cou- 
pling constant with time- like argument, leading to a more 
convergent series [9]. 

The last unknown scale in Eq. (14) can be estimated. 
It turns out that Qi ^ Q which is the value we have used 
[2] . The expressions for the coefficients 7^ and Hi can be 
found in Ref. [8], and the four- loops /3-function is given 
in Ref. [10]. The final result in numerical form in terms 
of a = cts/TT for QCD with five active flavors reads: 

^^e+e-^h^g) =1 + a{Qi) + 1.84a(Q2)' 

- 1.00a(Q3)^ - 11.31a(Q4)^ ■ (15) 

This is a more convergent result compared to previous 
estimates, and it is free of any scheme and scale ambigu- 
ities (up to strongly suppressed residual ones). 

To find numerical values for the effective scales, the 
asymptotic scale, A, of the running coupling must 
first be determined by matching Eq.(15) with exper- 
imental results [11]: jfRltp'^^iVs = 31.6 GeV) = 
1.0527 ± 0.0050 . Using a logarithmic expansion solu- 
tion of the renormalization group equation for a wc find: 
= 419tig| MeV. We have used the MS definition 
for the asymptotic scale, and the asymptotic scale of TZs 
can be taken to be the same for any S. The effective scales 
are found to be: Qi = 1.3 Q ,Q2 = 1.2 Q, Qa ~ 5.3 Q. 
The values are independent of the initial renormalization 
scale up to some residual dependence coming from the 
truncated /3-function, which is less than the quoted ac- 
curacy on the numbers. This is illustrated in Fig. 1. For 
we have taken the LO value, which is sufficient to 
get the conformal series at four loops. Its higher order 



value has artificial strong residual renormalization-scale 
dependence due to the large numerical value of 114 in 
QCD with five active flavors. These flnal scales deter- 
mine the effective number of quarks flavors at each order 
of perturbation theory [12]. 

For completeness, we use our final result to predict the 
strong coupling constant at the Z-boson mass-scale: 

a,(Afz)=0.132t°:^i? . (16) 

We can apply our result to Abelian QED, where 
j^e e -fh seen as the imaginary part of the QED 

four loop IPI vacuum polarization diagram by the opti- 
cal theorem, and find in this case nearly complete renor- 
malization scale independence of all three scales to the 
NNNLO due to the small value of the coupling constant. 
Numerically, we get for three (lepton) flavors: 

^^QED^'(O) =1 + 0.24a,(Qi) ~ 0.08ae(Q2)' 

- 0.13ae(Q3)' + 0.05ae(g4)* , (17) 

where = e^/Air and {^, ^, ^} = {1.1,0.6,0.5}. 

In this letter we have shown that a generalization of 
the conventional MS-scheme is illuminating. It enables 
one to determine the general (and degenerate) pattern 
of nonconformal {/3i}-terms and to systematically deter- 
mine the argument of the running coupling order by order 
in pQCD, in a way which is readily automatized. The re- 
summed series matches the conformal series, in which no 
factorially divergent n!/?"a" "renormalon" series appear 
and which is free of any scheme and scale ambiguities. 
This is the final expression one should use for physical 
predictions. The method can be applied to processes 
with multiple physical scales and is consistent with QED 
scale setting in the limit Nc — >■ 0. The new method sat- 
isfies all of the principles of the renormalization group, 
including the principle of maximum conformality, and it 
eliminates an unnecessary source of systematic error. 
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FIG. 1: The final result for i?" " as a function of 
the initial renormaUzation scale (solid line), 
demonstrating the initial scale-invariance of the final 
prediction, up to strongly suppressed residual dependence. 
The shaded region is the experimental bounds with the 
central value given by the dashed line. 
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